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Remarks on the chemical composition of highest-energy cosmic rays 
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Abstract: We present arguments aiming to reconcile the apparently contradictory results concerning the chemical com- 
position of cosmic rays of highest energy, coming recently from Auger and HiRes collaborations. In particular, we argue 
that the energy dependence of the mean value and root mean square fluctuation of shower maxima distributions observed 
by the Auger experiment are not necessarily caused by the change of nuclear composition of primary cosmic rays. 
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1 Introduction 

The identities of highest-energy cosmic rays remains still 
an open question. Possible conclusions on either protons 
or iron nuclei dominance in cosmic ray flux leads to prob- 
lems [1]. Seeking to determine the nuclear identities of 
ultrahigh-energy cosmic-ray (UHECR) particles, the devel- 
opment of extensive air showers (EAS) of secondary par- 
ticles in the atmosphere was extensively examined. The 
Auger collaboration [2] has determined both the shower 
maximum (X max ) (the penetration depth in the atmo- 
sphere at which the shower reaches its maximum number 
of secondary particles) and the complementary observable 
a(X max ) (the root mean square fluctuation of X max from 
event to event). Their results seem to indicate a transition, 
at primary energies of a few times 10 18 eV, from the flux 
dominated by protons to the one increasingly dominated at 
higher energies by iron nuclei. The HiRes collaboration [3] 
has analyzed event-by-event fluctuations of data in terms of 
the truncated fluctuation widths ar ( the X max distribution 
was truncated at 2a(X max ) ), and reaches a different con- 
clusion. We would like to present here arguments that the 
observed energy dependence of (X max ) and a{X max ) by 
Auger experiment are not originated by the changes of nu- 
clear composition of primary cosmic rays (cf., also, [4]) 
and that the highest-energy cosmic rays seems to be domi- 
nated by protons. 

2 Inconsistency in the iron abundance 

With the energy increase, the spectacular Auger data [2] 
show almost monotonic changes from proton composition 
towards iron one for both (X max ) and cr{X max ) observ- 
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Figure 1 : The energy dependence of relative abundance of 
iron in CR as extracted from (X max ) and a(X max ) given 
by Auger experiment [2] (in frame of QGSJETII [10] and 
EPOSvl.99 [11] models). 

ables. For (X max ) such dependence can be easily inter- 
preted by two component cosmic ray composition (with 
relative abundance of iron nuclei a and contribution of pro- 
tons 1 — a ), for which we expect 
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where (X max ) P and (X ma x)Fe are the shower maxima 
for pure proton and iron nuclei, respectively. However for 
a(X max ) we have nonmonotonic dependence on a, 

<r 2 = (1 — a) Cp + aa 2 Fe + 
+a (1 — a) ((X 

max) p {X max ) . (2) 

For this reason the experimental data (with similar en- 
ergy behavior) lead to quite different chemical composi- 
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tion, ranging from the proton dominated for (X max ) to the 
iron dominated for a(X max ) (cf. Fig.l). 

3 Importance of the first interaction point 

Some remarks are in order at this point (cf. [5, 6, 7, 9]). 
Most of the charged particles in the shower are electrons 
and positrons coming from the electromagnetic subshow- 
ers initiated by photons from ir° -decay, with energies near 
the critical energy (e = 81 MeV in air). The mean depth 
of maximum for an electromagnetic shower initiated by a 
photon with energy E 1 is 

{X™ x {E 1 ))=X \n{E 1 /e), (3) 

where Xo ~37 g/cm 2 is the radiation length in air. A 
nuclear-initiated shower consists of a hadronic core feed- 
ing the electromagnetic component primarily through tt° 
production. In general, for an incident nucleus of mass A 
and total energy E (including protons with A=\) the depth 
of maximum is expressed by 

(X max (E)) = {X™ x {{E/ A) {K/ („)))) + (Xi) , (4) 

where (X\) is the mean depth of the interaction with maxi- 
mal energy deposition into shower (usually called the depth 
of the first interaction), K denote inelasticity and (n) is re- 
lated to the multiplicity of secondaries in the high-energy 
hadronic interactions in the cascade. If the composition 
changes with energy, then {A) depends on energy and 
(X max ) changes accordingly. The situation is, however, 
essentially more complicated. Whereas for a primary nu- 
cleus in which the energy is to a good approximation sim- 
ply divided into A equal parts, in a hadronic cascade there 
is instead a hierarchy of energies of secondary particles in 
each interaction, and a similar (approximately geometric) 
hierarchy of interaction energies in the cascade. In this case 
(n) has to be understood as some kind of "effective" multi- 
plicity, which does not have a straightforward definition in 
general. For this reason the change of primary composition 
or the violation of Feynman scaling are widely discussed 
since many years. In addition to this, the inelasticity (K) 
can itself be function of energy [15]. 

The probability of having the first interaction point of a 
shower, X\ , at a depth greater than X is 

P(X 1 > A) ~ cxp {-X/X) , (5) 

where A is the interaction length. In the case of perfect cor- 
relation between X max and X\, i.e., when fluctuations in 
the shower development were nonexistent, one could use 
directly the exponential distribution of showers with large 
X m ax to calculate X\ and hence the proton-air cross sec- 
tion. However, intrinsic shower fluctuations modify rela- 
tion between the depth of maximum distribution and the 
interaction length. This modification is typically expressed 
by a factor k = A/A and leads to P(X max > X) <~ 
exp(— X/X). The factor k depends mainly on how fast 



is the energy dissipation in the early stages of shower evo- 
lution. In particular it is sensitive to the mean inelasticity 
and to its fluctuations. In general, a model with small fluc- 
tuations in secondary particle multiplicity and inelasticity 
is characterized by a smaller k factor than a model with 
large fluctuations. Under the assumption of similar fluc- 
tuations in multiplicity and inelasticity, a model predict- 
ing a large average number of secondary particles leads to 
smaller overall fluctuations of the cumulative shower pro- 
file of the secondary particles and hence to a smaller k fac- 
tor. 

In the absence of internal fluctuations, all showers would 
develop through the same amount of matter, AX = 
X m ax — Ai, between the first interaction point and the 
maximum. As a consequence, a perfect correlation be- 
tween Xmax and Ai would exist, and their distributions 
would have exactly the same shape, shifted by a constant 
AA. In that case the slope of the X ma x distribution, A, 
would be equal to the mean interaction length, A. Intrinsic 
fluctuations in shower development (after the first interac- 
tion) affect the relation between the interaction length A 
and the slope A that describes the exponential tail of the 
Xmax distribution. The relation is often expressed with a k 
factor k = A/ A. For more properties of EAS and influence 
of shower fluctuations on studies of the shower longitudinal 
development see Ref. [5, 6, 7, 8, 9]. 

The effect of fluctuations in AA is to broaden the corre- 
lation of Xmax with Ai . However, we can roughly write 
that 

a(Xmax)=<j(X 1 )+t(<T(AX)), (6) 

where a (Ai) oc< Ai > and the function £ describes influ- 
ence of shower fluctuations after the first (main) interaction 
point (notice that for the probability distribution given by 
Eq.( 5) the fluctuation in Ai is er(Ai) = \JVar{X\) = 
(Ai), whereas for Ai interpreted as the main interaction 
point we have er(Ai) = (Ai) / ' y/R where k determines 
in which of the succesive interactions of projectile particle 
the energy deposition to the shower is maximal). Because 
of Eq.(4), where (X ma x) = {X™ x ) + (*i>. we can con- 
struct observable in which influence of fluctuation of the 
first interaction point is strongly suppressed, namely 

{Xmax) ^{Xmax) = 

= {XZZx {{E/A){K/ < n >)> + £ (er(AA)) .(7) 

4 Results 

In Fig. 2 this observable is plotted for Auger [2] and HiRes 
[3] data in comparison with different models [10, 11, 12, 
13]. To make the results from both experiment to coincide, 
the HiRes data are shifted by 10 g/cm 2 (in this case pre- 
dictions from QGSJETII model are roughly the same for 
both experiment) 1 . 

1. Notice that Auger compares the data with pure simulations. 
HiRes quotes data including all detectors effect and compare it 
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Figure 2: (X max ) — cr(X max ) as deduced from Auger data and (X max ) — aT{X max ) from HiRes data in comparison 
with different models. Notice that HiRes experimental data are shifted by 10 g/cm 2 to make the model predictions for 
both experiment coincide. 



Notice that (X max ) — cr{X max ) is still dependent on mod- 
els and, in particular, it is sensitive to the chemical compo- 
sition. Showers initiated by protons are seemingly different 
from those initiated by iron nuclei. 

From Fig. 2 we can learn that the chemical composition is 
not the origin of the effect observed by Auger experiment. 
Moreover, the experimental data fairly well coincide with 
the proton dominant primary composition. Within the toy 
model of primary composition (only two components: iron 
nuclei with relative abundance a and protons with abun- 
dance 1 — a ) we can evaluate a from (X max ) — a(X max ) 
as given by Auger experiment. The results is shown in 
Fig. 3. For the reference model QGSJETII the abundance 
of iron is roughly independent on energy (a ~ 0.05 -j- 0.1) 
and even for model EPOS v. 1.99 [11], which leades to the 
maximal abundance of iron, it increases slowly with energy 
(varying in interval a ~ 0.15 4- 0.3). The iron abundance 
shown in Fig. 3 coincides with the one which can be esti- 
mated from HiRes data. The comparison of a from Auger 
and HiRes data is shown in Fig. 4. In the energy region 
2 • 10 18 ~ 5 • 10 19 eV the mean values of a, evaluated 
from (X max ) — cr{X max ), are equal a = 0.08 ± 0.01 from 
Auger data and a = 0.06 ± 0.05 from HiRes data (no- 
tice that HiRes data on (X max ) result in comparable value, 
a = 0.03 ±0.02). 

5 Possible interpretation 

From Fig. 2 we can learn that (Xi(E)) gives the main 
contribution to the energy dependence of (X max ) and 
a(X max ) observed experimentally. Two factors can af- 
fected energy dependence of (Xi(E)): the cross sec- 
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Figure 3: The energy dependence of relative abundance of 
iron in CR as extracted from (X max ) — a(X max ) as given 
by Auger experiment and shown in Fig. 2. 

tion (interaction mean fee path A) and the inelasticity K. 
Roughly, (Xi) = X ■ k, where n determines in which of 
the successive interactions of projectile the energy depo- 
sition to the shower is maximal. For a uniform inelastic- 
ity distribution in the maximal possible interval for a given 
(A") one has k ~ 1 + 1.85(0.75 - (K)). The rapid in- 
crease of inelastic cross section in energies E > 10 18 
eV cannot be excluded. In particular, if gluon saturation 
occurs in the nuclear surface region, the total cross sec- 
tion of proton— nucleus collisions increases more rapidly 
as a function of the incident energy compared to that of a 
Glauber-type estimate [14]. Although in [15] the decrease 

to the models 'after' the detector simulation. Unfortunately that 
means that both approaches cannot be compared directly. 
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Figure 4: The energy dependence of relative abundance of iron in CR as seen from Auger and HiRes data (in the frame 
ofQGSJETII model [10]). 



of inelasticity (K) with energies was discussed in lower en- 
ergy region, its increase at energies E ~ 10 18 eV is be no 
means excluded (cf. the percolation effects which in high 
energies leads to increase of inelasticity [16]). Both possi- 
bilities are questionable and require an abrupt onset of new 
physics beyond the standard model (notice however that 
here, the center of mass collision energy is about few hun- 
dreds of TeV, far beyond that can be studied at LHC). Tak- 
ing into account the HiRes data (where X max distribution 
was truncated at 2a ) we can learn that the tails of X max 
distribution are crucial. For this reason, the role of biases 
due to the small statistics in analyzing CR data of highest 
energy remains an open question (cf. ref. [4]). It is inter- 
esting to note that the observable (X max ) — <r(X max ) is 
rather insensitive to the possible biases of the tail of X max 
distribution [4]. 

6 Concluding remarks 

To summarize, we conclude that the spectacular energy de- 
pendence of the shower maxima distribution reported by 
Auger collaboration [2] is not necessarily (or not only) due 
to the changes of chemical composition of primary cos- 
mic rays. The observed effect seems rather to be caused 
by the unexpected changes of the depth of first interac- 
tion in energies above 2 10 18 eV. They would requires, 
however, an abrupt onset of some "new physics" in this 
energy region and are therefore questionable. We argue 
that it would be highly desirable to analyze the observable 
(X max ) — cr(X max ) in which fluctuations of the depth of 
the first interaction, as well as the possible biases of the 
tail of X max distribution, are strongly suppressed. This 
observable still depends on the model of multiparticle pro- 



duction and is sensitive to the chemical composition of the 
primary cosmic rays. 
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